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ABSTRACT 

Finite source effects can be important in observations of gravitational microlensing of stars. Near caustic 
crossings, for example, some parts of the source star will be more highly magnified than other parts. The 
spectrum of the star is then no longer the same as when it is unmagnified, and measurements of the atmospheric 
parameters and abundances will be affected. The accuracy of abundances measured from spectra taken during 
microlensing events has become important recently because of the use of highly magnified dwarf stars to probe 
abundance ratios and the abundance distribution in the Galactic bulge. While the abundance ratios in general 
agree with the giants, with the possible exception of [Na/Fe], there may be more dwarfs than giants with 
high [Fe/H]. In this paper, we investigate the effect of finite source effects on spectra by using magnification 
profiles motivated by two events to synthesize spectra for dwarfs between 5000K to 6200K at solar metallicity. 
We adopt the usual techniques for analyzing the microlensed dwarfs, namely, spectroscopic determination of 
temperature, gravity, and microturbulent velocity, relying on equivalent widths. We find that ignoring the finite 
source effects for the more extreme case results in errors in T e ff< 45K, in log g of <0.1 dex and in £ of <0.1 
km/s. In total, changes in equivalent widths lead to small changes in atmospheric parameters and changes 
in abundances of <0.06 dex, with changes in [Fe I/H] of <0.03 dex. For the case with a larger source-lens 
separation, the error in [Fe I/H] is < 0.01 dex. This latter case represents the maximum effect seen in events 
whose lightcurves are consistent with a point-source lens, which includes the majority of microlensed bulge 
dwarfs published so far. 

Subject headings: Galaxy: abundances, bulge - gravitational lensing - stars: atmospheres - techniques: spec- 
troscopic 
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1. INTRODUCTION 

A microlensing event occurs when a stellar mass object 
passes between the observer and a background source. For 
events in the Local Group, the source is a star. If the source 
is finite, rather than infinitely small, parts of the source may 
be magnified more than other parts. In the case of stars, for 
example, the limb can be brightened relative to an unmagni- 
fied source. Since the light from the limb comes from dif- 
ferent temperature profiles in the star c ompared to the cen- 
ter, this will affect the spectrum (e.g. IValls-Gabaudl 1 1 9951 
lLoeb & Sasselovl[l995l) For example, a spectrum of the limb 
of the Sun has Balmer lines with weaker wings, while lines of 
other elements ca n be either weakened or strengthene d rela- 
tive to the center (lHastingslll873HHale & Adamslfl907l) . 

Recently, observers using large telescopes have published 
abundance ratios in dwarfs in the Galactic bulge from spectra 
obtained while the sources were highly magnified. A num- 
ber of these have turned out to have super- solar metallici - 
ties, including OGLE-20 06-BLG-265S (IJohnson et aTJl2007l) . 
MOA-2006- 0BLG-099S (IJohnson et al.ll2008l) . OGLE-2007- 
BLG-349S (ICohen et all 2008). M OA-200 8-BLG-310S and 
MOA-2008- BLG-3 1 1 S (ICohen et alJl2009l) and OGLE-2007- 
BLG-514S dEpstein et al.1 120091) . Some are metal-poor in- 
cluding the subgiant OGLE-2008-BLG-209S (Bens bv et ail 



2009a) and the dwarf OGLE-2009-BLG-076S dBensby et al.1 
2009b), which, at [Fe/H]= -0.76 is the most metal-poor 



dwarf/subgiant observed when microlen s ed. Fo r a sample of 
eight microlensed dwarfs, lEpstein et al.1 (120091) found that a 

1 Current address: Institute for Advanced Study, 1 Einstein Drive, Prince- 
ton, NJ 08540, USA 

Electronic address: jaj.dong, gould@astronomy.ohio-state.edu 



K-S test gave a 1.6% chance that the dwarfs were drawn from 
the same metallicity distribution function (MDF) as th e giant 
MDF from Zoccali et al. (2008). Bensby et al. (2009c) report 
results for a total of 13 dwarfs, including new results for more 
metal-poor dwarfs, so clearly the compar ison between the gi- 
ants and the dwarfs is an evolving topic. ICohen et al.l (2008) 
proposed that mass-loss on the giant branch prevents some 
more metal-rich stars from b ecoming red giants , similar to 
the mechanism suggested by Kali rai et al.1 (120071) to explain 
the low-mass He white dwarfs in the metal-rich cluster NGC 
679 1 . However, this explanation predicts that the red clump 
MDF is biased to more metal-poor stars because of metal- 
rich stars skipping the horizontal branch phase and that the 
red giant luminosity function should drop more than expected 
from theoretica l models, predictio ns that do not agree with 
the giant data dZoccali et al.1 12008). Another explanation is 
that there are system atic differences bet ween giant and dwarf 
abundance analyses. ICohen et al.1 (12009) calculated that a sys- 
tematic offset of 0. 10 dex between the dwarf and giant metal- 
licity scales would give the mean metallicity offset a signif- 
icance of 2-cr and larger systematic offsets wo uld obviously 
decrea se the significance further still. Finally, IZoccali et al.1 
(2008) suggested that the spectra of microlensed dwarfs could 
be affected by differential magnification sufficiently that the 
usual analysis of the dwarfs, which does not take this into 
account, could lead to biased answers and help explain the 
possible discrepancy. This last suggestion can be tested by 
comparing the answers obtained from synthetic spectra with 
and without differential limb magnification. 

While most microlensing events follo w the lightcurve of 
a point source, about ^3% dWittl 1 1995D of events show fi- 
nite source effects. This fraction is even higher for high- 
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magnification events that are targets of the current generation 
of dwarf studies. Of the eight published events for which 
dwarf spectra were obtained, we know that at least three 
of the events were affected by finite source effects, namely 
OGLE-2007-BLG-514, OGLE-2007-BLG-349, and MOA- 
2008-BLG-310. We need to determine the size of the effect 
that differential limb magnification (DLM) has on the spec- 
tra and the measurement of the effective temperature (T e ff), 
gravity (log g), metallicity ([Fe/H]), microturbulent velocity 
(£) and abundance ratios to correctly interpret these events. 

In addition to probing the chemical evolution of the bulge, 
the accuracy of the measured T e ff from the spectrum is impor- 
tant for testing the method by which colors of source stars in 
microlensing events are determined. Using a metallicity and 
a T e ff derived from the spectrum, we can predict the color of 
th e star using relations between color and T e f t , such as that 
by iRar mrez & Melendez (2005), and compare with the color 
estimated using standard microlensing techniques, which rely 
on the offset of the star from the red clump. The results so far 
indicate that if the color of the red clump is (V-/)o=1.05, the 
T e ffS derived spectra are in agreement with T e ffS from colors. 

Much work has been done on the effects of DLM of the disk 
of giants during a microlensing event because these events 
are easier to find and are longe r-lasting than similar events 
in dwarfs. Vall s-Gabaudl (1 19981) predicted the effects on the 
spectrum of a giant, in particular the H/3 and CO lines, by 
fitting simple analytic expressions to the limb p rofiles for 
both continuum and lines from the iKuruczl (1 19921) models. 
He found that the equivalent widths could change by 20% 
over the course of the event. iHeyrovsky et al.l (120001) com- 
puted more realistic models avoiding the use of linear approx- 
imations. They calculated contribution functions for lines, 
which indicate where the lines are formed, and illustrated 
the range of behavior individual spectral fe atures will have 
becaus e of varying center-to-limb profiles. iGaudi & Gould! 
( 1999) pointed out that binary lens events had particular ad- 
vantages for resolving stellar surfaces, for example the fact 
that the second caustic crossing can be predicted if the event 
is monitored, and therefore carefully observed. 

Thanks to intensive monitoring by observers, the effect of 
the size of the source on microlensing events has been ob- 
served many times. Finit e source effects we re observed for 
the giant MACHO 95-30 (lAlcock et al.ll 19971) . including vari- 
ations in the equivalent widths of Ha and the TiO bands. The 
bulge K3 giant star EROS BLG-2000-5S was highly magni- 
fied by a binary lens. Because the caustics of a binary lens 
are symmetric, the timing of the second caustic crossing was 
predicted to high accura cy and allowed int ensi ve follow-up 
during this event. Both lAlbrow et all (1200 ll) and lCastro et all 
(2001) found decreases in the strength of the Ha line when 
the limb was more highly magnified relative to the center. 
lAfonso eTa l. (2001) monitored the same event photometri- 
cally, measured the color of the limb, and found that the 
outer 4% of the limb showed strong Ha emission, which they 
ascribe d to the presence of a chromosphere. Albr ow et al.l 
( 1999) performed the first limb-darkening measurement of a 
bulge giant for the event MACHO 1997-BLG-28, using pho- 
tometry in the V and / bands. Other limb-darkening mea- 
surement s were done for gia nts including MACHO 1997- 
BLG-41 (lAlbrow et all 12000), and OGLE-2002-BLG-069S 
(iKubas etal.1 120051) . which the authors found agreed with a 
linear limb-darkeni ng law. An opposite result was found by 
ICassan et al.l d2006) who obtained photometry for the K gi- 
ant OGLE-2004-BLG-254S and derived limb darkening co- 



efficients in the / and R bands. Combining these results 
with other microlensed-based limb darkening measurements 
of K giants and comparing them to the predictions of AT- 
LAS9 models, they found a disagreement, which they sug- 
gested was the result of the inadequacy of the linear limb 
darkening law. Time -resolved spectra o f OGLE-2002-BLG- 
069S, a G5III star, bv lCassan et al.l d2004l) . showed the power 
of this technique for probing the stellar atmosphere. In this 
case, the changes during the event in wings of Ha line, which 
are formed deep in the atmosphere, agreed with model pre- 
dictions while the core of the line showed an emission peak, 
which again can only be explained by a chromosphere. That 
the temperature structure in the outer layers is not well un- 
derstood for this event was also found by the analysis of 
iThurl et ail d2006l) . These studies show the potential to in- 
vestigate the atmospheric structure of giants. 

Stud ies of dwarf st a rs hav e been much rarer. Afon so et al.l 
(2000) and I Abe et al.l (120031) reported the only measurements 
of limb darkening in dwarfs, for MACHO 98 SMC-1, a main- 
sequence A star and MOA 2002-BLG-33S, a solar-type star, 
respec tively. The most relev ant work to this paper is the anal- 
ysis of iLennon et all (fl996l) of the event MACHO 96-BLG-3. 
They took three exposures when the source, a dwarf star, was 
moving across a caustic. These R^llOO spectra with signal- 
to-noise (S/N) of 25- 100 did not show an y convincing cases 
of profile variability. ILennon et al.l (1 19961) compared the ob- 
servations of this star with both a library of high-resolution, 
high S/N spectra of F and G stars observed at Calar Alto and 
a grid of synthetic spectra, and derived stellar parameters of 
T eff =6 100K from the Ha line, log §=4.25 from the Mg I triplet 
and a metallicity ([M/H]) between 0.3 and 0.6 from fits to 
regions with many Fe I and Ca I lines. They calculated the 
expected deviations from the unmagnified spectrum for this 
event, and found that the expected change in the line profiles 
for the three spectra was < 1%, while the expected change in 
the continuum was < 2%. Given the small changes expected 
and the S/N of the spectra, it is not surprising that no changes 
were observed, and that the derived atmospheric parameters 
of the star would also not be affected. Indeed, it is not surpris- 
ing that most changes are small. The amount of magnification 
depends on the distance from the lens, but the emitted spec- 
trum is the same for an entire annulus (for a spherical star). 
Because the distance from the lens varies around the annulus, 
the average magnification of a spectrum at a particular an- 
nulus is smaller than the largest magnification of an particular 
spot would suggest. In addition, while the spectrum of the star 
increasingly changes from center to limb, the intensity of the 
limb is lower than the center, by factors of a few. Therefore it 
is difficult to overcome the influence of the light coming from 
the central regions of the star in the disk-averaged light. 

The work discussed above focused on learning about the 
atmosphere of the star and showed that finite source effects 
cause changes in the spectrum of the star, which can be in- 
verted to give temperature as a function of depth. Because 
the same effects that make these extreme events interesting 
as probes of the atmosphere will also change the spectrum 
that is analyzed for elemental abundances, the goal of this pa- 
per is to determine the effect of DLM on the T e ff, £, log g, 
[Fe/H] and abundance ratios using the standard spectroscopic 
techniques, such as equivalent width analysis. Therefore, we 
will be concerned with mimicking standard abundance analy- 
sis as closely as possible, rather than exploring the full range 
of knowledge that can be gleaned from time-resolved spec- 
troscopy of microlensed events. 
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FIG . 1 . — Differential limb magnification profiles for Event A. The bold line 
shows Profile 1 , which represent the most extreme differential magnification 
for this event, and is the profile from the first 15 minutes of Event A. At most, 
the limb is magnified by ~ 30% more than the center. 

2. MAGNIFICATION PROFILES DURING 
MICROLENSING EVENTS 

Our work was motivated by two events, MOA-2008-BLG- 
310 and OGLE-2007-BLG-514, which are illustrations of the 
kind of effects seen. For MOA-2008-BLG-310, the lens did 
in fact transit the source, and there were small perturbations 
durin g this transit caused by a planetary companion in the 
le ns jj anczak et al. 2009). However, the spectrum analyzed 
bv lCohen et al I d2009l) began at UT 22:51, which was 21 min- 
utes after the end of the transit, when the lens and source cen- 
ter were separated by z = 1.22 source radii. Figure Q] shows 
the annulus-averaged profiles for a point lens at 15 minute in- 
tervals for the MOA-2008-BLG-310 geometry, beginning at 
z=1.2 and continuing as the lens and source moved further 
apart. The magnification is given as a function of 9, which is 
the angle between the normal to the stellar surface and the line 
of sight to the observer. To determine This look at the later 
stages of MOA-2008-BLG-310 will be referred to as Event 
A. The start time for t he first profile is three minutes before 
the lCohen et al.l d2009l) observations began, and thus serve as 
a direct measurement of the size of the effect of DLM on those 
abundances, but they also serve a broader purpose. 

In this particular case, we know that z = 1.22 at the time 
of the observations because the source size was earlier de- 
tected by observations (from Africa) of a direct source cross- 
ing. However, if the closest separation between source and 
lens during the entire event (i.e., the impact parameter) had 
been zo = 1.22, there would have been no source crossing at 
any time. If the impact parameter is sufficiently large, then it 
becomes impossible to determine the source radius from the 
lightcurve. In such a case, only upper limits can be placed on 
the magnitude of DLM and profiles of the DLM as a function 
of 8 cannot be derived. 

However, for zo < 1 .2, the effects of the finite-source size 
on the light curve are sufficiently pronounced to measure zo- 
Thus, the top curve in Figure Q] represents the most extreme 
case of DLM that would occur without being noticed, and 
thus serves as a general check on ignoring DLM when it is 
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FIG. 2. — The top panel shows, for Event B, the source-averaged magnifi- 
cation as a function of time, measured in units of the source-radius crossing 
time. The bottom panel shows the geometry of the event at a moment when 
the source (circle) is just exiting the caustic (acute-angled structure), which 
is a contour of formally infinite magnification. The source-center position is 
marked by small circles at 5 minute intervals. As the source passes over the 
caustic, it is differentially magnified, with first the limb, then the center, and 
finally the limb being the most magnified. 

not detectable from the light curve. 

The other case represents a more extreme event, inspired by 
OGLE-2007-BLG-514, and will be called Event B. Here, the 
source trajectory crossed a cusp from a binary lens and pro- 
duced extreme magnification variations (Figure O. Once the 
parameters are selected, it is straightforward to compute the 
magnification at each time and for every point on the source 
plane. We calculated the mean magnification in concentric 
rings with annuli eq ual to 0.01 sour ce radii using the "loop- 
linking" technique dDong e t al. 2006). The profiles are spaced 
at five minute intervals for this event. Actual observations are 
always longer, usually 5-6 x longer in order to get enough S/N 
in each exposure to reliably extract the spectrum. Spectra at 
5 minute intervals could only be obtained for dwarfs mag- 
nified to apparent magnitudes that have not been seen in an 
event to date, so in reality microlensing spectra will smear out 
these profiles and dilute their effects. Figure [3] shows 51 dif- 
ferent profiles that occurred during this event. We note that 
in this case we have ample notice from the lightcurve that fi- 
nite source effects are important and we could use that infor- 
mation to derive the magnification profile at the time spectra 
were taken and make model spectra that include the effects 
of DLM. However, in this paper, we consider the cases that 
ignored finite source effects on the spectrum to measure the 
size of errors that this induces on the parameters, metallicities 
and abundance ratios. 

3. SYNTHESIZED SPECTRUM 

We explored the changes to the spectra of dwarfs caused 
by the magnification profiles given in §2. We interpolated 
a set of model atmospheres from the K urucz-Castelli grid 
with n ew opacity distribution functions (Castel h~& Kuruczl 
d2003l) 2 . The models are solar- me tallicity and are spaced 

2 Available at http://wwwuser.oat.ts.astro.it/castelli/grids.html 
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FIG . 3 . — The magnification of the disk of a dwarf as a function of radius at 
each interval marked in Figureff] We see in the most extreme cases that parts 
of the disk are magnified up to 5 X more than other parts. The two profiles 
that are examined in more detail in the text are for time step 18 (bold line 
peaking at 9 = 0) and time step 35 (bold line peaking at 8 = 0.75) 



TABLE 1 
Parameters for Model 
Atmospheres 



Teff(K) 


logs 


[Fe/H] 


£ (km/s) 


5000 


4.60 


0.00 


1.500 


5100 


4.59 


0.00 


1.500 


5200 


4.57 


0.00 


1.500 


5300 


4.55 


0.00 


1.500 


5400 


4.53 


0.00 


1.500 


5500 


4.51 


0.00 


1.500 


5600 


4.48 


0.00 


1.500 


5700 


4.46 


0.00 


1.500 


5800 


4.43 


0.00 


1.500 


5900 


4.39 


0.00 


1.500 


6000 


4.33 


0.00 


1.500 


6000 


4.06 


0.00 


1.500 


6100 


4.01 


0.00 


1.500 


6200 


3.98 


0.00 


1.500 



every 100 K from 5000K to 6200K. This temperature range 
covers the range of temperatures for dwarfs observed in the 
bulge. Cooler dwarfs (< 5000 K) are unlikely to ever be tar- 
geted because of their intrinsic faintness and because if the 
dwarfs are as cool as the giants, many of the advantages of 
measuring abundances in dwarfs rather than giants (such as 
reduced blending) no longer apply. Although the majority of 
the bulge population is old (lOrtolani et al.ll 19951) . we wanted 
to explore a large temperature range of main-sequence/main- 
sequence turnoff stars, and th erefore we ado pted log g values 
from a Yale- Yonsei isochrone ( Yi et al.l20 1 ) of 4 Gyr. At this 
age, stars with T e ff=6000K have two possible values for log 
g, because they are at the turnoff, while dwarfs with 6100K 
and 6200K are still present because of the blue hook in the 
isochrone. The two 6000K dwarfs allow us to measure the ef- 
fect that small changes in log ghave on the resulting spectra. 
Table 1 lists the temperatures and log g's for the dwarfs. 



3.1. Method of Synthesizing Spectra 

We focused our attention on three 200 A sections of the 
spectrum, centered on Ha (6460 A- 6660A), H/3(4757A- 
4957A) and the Mg I triplet (5067 A- 5267 A). These regions 
have both strong lines, which are expected to show the largest 
variations, as well as a number of weaker lines whose equiv- 
alent widths would be used in an abundance analysis. The 
wavelength range from blue to red also ensures that detectable 
unblended lines can be found for both cooler stars, which have 
crowded blue regions, and hotter stars, which tend to have 
weak lines in the red. The Ha and H(3 lines are also used as 
a temperature indicator in hotter d warfs. For each model at- 
mosphere, we used Turbospectrum dAlvarez & Pled [l998) to 
generate intensities 1(6') at 100 values of 9. To determine the 
total flux from the star, we added up the intensities coming 
from annuli from the center to the radius R of the star 



F = 2n 



I(0)rdr 



(1) 



(see e.g jMihalasll 19781) Using r = RsinO, for the unmagnified 
case, the observed flux can be obtained by numerical integra- 
tion of the equation 



F = 2tt 



I(9)cos9sin9d6 



(2) 



This gives the same answer as when Turbospectrum outputs 
a flux, rather than intensities. For the magnified cases, each 
annulus was multiplied by its magnification factor before in- 
tegration. The spectra were smoothed to a FWHM of 0.1 1 A, 
or R~ 45,000, similar to the resolution at which the dwarfs 
are observed. 

3.2. Comparison of Spectra 

To investigate the size of possible systematic effects on the 
abundances in dwarf stars, we wished to determine the max- 
imum effect on the spectroscopic analysis, and therefore be- 
gin by identifying the cases for which the spectra deviate the 
most from the unmagnified case. We took the ratio between 
the magnified and the unmagnified spectra, renormalized the 
spectra and then calculated the rms. 

For Event A, neither the profiles nor the rms varies much, 
but the largest deviations are found for the profile calculated 
for the first time step, indicated by the bold line in Figure [T] 
For Event B, the rms is highest, as expected, when the limb is 
magnified by a high factor or when the center region is mag- 
nified the most (bold lines in Figure [3). These occur for the 
1 8th and 35th profiles calculated, corresponding to 85 and 1 70 
minutes after Event B began. We will use these three profiles 
as examples to calculate the size of the effects on the abun- 
dances. The indicated case from Event A will be called Pro- 
file 1 and the two indicated cases from Event B will be called 
Profile 18 and 35. 

Figure |4] illustrates the changes to a spectrum by showing 
the ratio of the unmagnified spectrum to the spectra created 
using the three different magnification profiles for the 5500K 
dwarf. The changes are small, but vary from line-to-line based 
on the location of the line-forming region. In general, the 
strong lines, including Ha and H/3, become weaker for Pro- 
files 35 and 1, for which the limb is magnified relative to the 
center, while the reverse is true for the spectrum from Profile 
18. The opposite behavior is seen for most of the weak lines. 
As expected, the deviations are largest for the two cases where 
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the magnification profile is more extreme and much smaller 
for Event A, where the differential magnification effects were 
~ 30%, rather than ~ 500%. 

4. ABUNDANCE ANALYSIS 

The goal of this paper is to determine the practical effect 
of finite source effects on abundance measurements in mi- 
crolensed dwarfs. To do this, we followed the usual steps 
used to analyze these spectra. Because differential redden- 
ing across the face of the bulge leads to uncertainty in the 
colors and magnitudes of stars in the bulge, analyses of mi- 
crolensed dwarfs have relied on spectroscopic methods of 
deriving atmospheric parameters, rather than color or appar- 
ent magnitude-based methods. DLM can impact the derived 
abundances in the source star by changing the EWs of the 
lines and by changing the atmospheric parameters derived 
from those EWs. To examine these effects and compute the 
total effect of DLM on the abundances, we first measured 
EWs in the synthetic spectra for the unmagnified case and 
Profiles 18, 35 and 1. We ran all sets of EWs through the 
original model atmosphere to see what changes. We found 
small changes in the diagnostics used to determine T e ff, log 
g, and £, in addition to changes in the abundances. We first 
discuss the magnitude of the changes demanded by the magni- 
fied spectra on the atmospheric parameters, considering each 
one in isolation. While the exact magnitude of the changes 
depends in detail on which lines are used for any particular 
study, the results here will give a general indication. Looking 
at the changes in each parameter separately shows the size 
of the effect if only some of the parameters are determined 
from the spectra. However, because deriving the atmospheric 
parameters usually depends on the other parameters, we next 
consider the total effect of the accumulated atmospheric pa- 
rameter changes plus EW changes. Once again, depending on 
the line list used, different correlations between the parame- 
ters may be found, but the overall tendencies can be general- 
ized. 

4. 1 . Equivalent Width Measurement 

In measuring EWs on the model spectra, we encountered 
the same concerns about continuum placement and blending 
as in measuring EWs on observed spectra. Because of our 
desire to determine differences caused by differential limb 
magnification, we focus on measuring EWs in a consistent 
manner. We first used IRAF 3 to do the continuum division, 
keeping the parameters of the fit, such as the order of the 
polynomial and the high and low-reject sigmas, the same for 
both magnified and unmagnified spectra. Then we measured 
the EWs without adjusting the continuum interactively, us- 
ing SPECTRE (C. Sneden, 2007, private communication) to 
fit Gaussian profiles to the s pectra. Our initial li nelist was se- 
lected using the solar atlas of Moor e et al.l(l 19661) as a guide to 
unblended lines, but this was not completely foolproof at the 
lower resolution (compared to the solar atlas) of the synthe- 
sized spectra, and over the large temperature range spanned 
by the models. Therefore we checked the list by determining 
the abundances from the lines for the unmagnified spectra. 
Any line that deviated by more than 0. 1 dex from the mean 
is eliminated. Also eliminated are lines with EW > 150mA 

3 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 



in the unmagnified case, because these lines are usually elim- 
inated in high-resolution analysis due to their sensitivity to 
damping parameters, continuum placement, blending in the 
wings and incorrect temperature structure in the outermost 
layers of the model atmospheres. For each temperature, the 
same lines were measured for each magnification profile; 
however, because the strength of lines varies with tempera- 
ture, the linelist changes for different temperatures. Because 
the same transition information, such as g/-value, was used 
both to synthesize the spectra and to calculate the abundance, 
there are no uncertainties introduced in this process by atomic 
data uncertainties. 

In Figure |3J we compare the difference in EWs for the three 
magnification profiles under study. We find that, as expected, 
the smaller differences between limb and center for Profile 
1 result in smaller changes in the EWs. As seen in Fig- 
ure 4, weaker lines are mostly strengthened and strong lines 
mostly weaken for Profile 35, while the opposite is true for 
Profile 18. This emphasizes the point that the magnitude of 
the changes for a particular analysis will depend on the dis- 
tribution of EWs that are being used. We then ran the EWs 
through Turbospectrum, using the original model atmosphere. 
This comparison isolates the effect of the EW changes on the 
abundances, before we consider changes caused by new at- 
mospheric parameters, and would reflect the total changes if 
the analysis relied on non-spectroscopic methods for atmo- 
spheric parameter determination, such as color of the unmag- 
nified star for T e ft and position on the CMD for gravity. The 
changes in abundances are shown in Figure [6] and show the 
magnitude of 5 loge is generally < 0.05 dex. 

4.2. r eff 

There are two standard ways of determining the tempera- 
ture of a star from its spectrum, without relying on its color. 
The first uses the wings of the Balmer lines, while the sec- 
ond requires excitation equilibrium for Fe I, so that there is 
no trend in a plot of lower excitation potential to abundance 
derived for each line. 

Ha is well-known to show the effects of DLM. Figure [7] 
compares the Ha profiles of the unmagnified case to the 
most extreme magnified cases (Profile 18 and 35) for the 
T e ff=5500K model. It is clear that the Ha profiles are different 
in the magnified models, and a different temperature would 
be derived if DLM were not taken into account. By compar- 
ing the magnified cases with the unmagnified profiles of the 
different temperature models, we determine that the T e g we 
would derive would be 75K or 100K hotter for Profile 18 and 
would be 75K or 100K cooler for Profile 35. This is consistent 
with the magnification profiles for these two cases, because 
Profile 18 has the center, where we see to hotter temperatures, 
magnified more than the limb, while the opposite is true for 
Profile 35. 

Determining the temperature via the wings of the Ha line 
requires high S/N data and becomes increasingly difficult at 
lower temperatures. Therefore, the more common method, 
used by all the papers on microlensed dwarfs discussed in 
the Introduction, is excitation equilibrium. We found that the 
slope of the line, determined by a least-squares fit, in the ex- 
citation potential-loge(Fe) plane changes if we used the EWs 
measured from the DLM spectra. Because of errors in con- 
tinuum division and blending, the slope of the line in the un- 
magnified case, which theoretically should be zero, deviates 
slightly from that value. Therefore, to calculate the change in 
T e ff from this method caused by the effects of DLM on the 
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series of model atmospheres with the temperature decreased 
by 25K, found the new slope using the magnified EWs, and 
used the results to calculate the necessary temperature change. 
The results are illustrated in Figure[8] We also did the test with 
+25 K and got the same answer to within 1 degree. 

4.3. Microturbulent Velocity 

The microturbulent velocity is determined for 1 -D analysis 
by ensuring that the slope for the derived abundance as a func- 
tion of reduced EW for lines of an element is zero. Usually, 
the large number of Fe I lines make it the element of choice 
for this test, and this was the case for the analyses of the mi- 
crolensed dwarfs. We compared the slopes of the lines in the 
Fe I vs. reduced EW (RW) plot in the magnified case to that 
in the unmagnified cases using the Fe abundances determined 
with the original model atmosphere in both cases. There is a 
change in the slope of the line in the magnified cases, because 
their EWs are affected by DLM.. As with the T e ff, the slope 
of the line in the unmagnified case is not exactly zero. We run 
our EWs and original model atmospheres through Turbospec- 
trum, this time with £ set to 1 .4 km/s and use the change in the 
slope for 0. 1 km/s to interpolate the change in £. The magni- 
tude of the change is always < 0.1 km/s (Figure |9). 



spectrum, we calculated the change in temperature needed to ° ° 

get the slope back to that measured in the unmagnified case. For the analysis of the microlensed dwarfs, gravity is deter- 

Because we found the temperature shifts are small, we ran a mined spectroscopically by demanding that Fe I and Fe II (or 
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other ions) give the same abundance. Obviously, changes in 
the EWs of these lines as well as changes in T e g and £ will 
result in a change in the derived log g. In Figure [TUJ we show 
the changes in log g that are necessary to get agreement be- 
tween Fe I and Fe II considering the changes in EW alone. 

4.5. Abundances 

In the previous subsections, we have considered the 
changes necessary in T e ff, £ and log g if we kept the other pa- 
rameters equal to those of the original model and only mod- 
ified the one under consideration. This is appropriate when 
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one model atmosphere parameter is chosen independently of 
the others. However, when deriving a model atmosphere, 
changes in one parameter will usually propagate and lead to 
total changes that can be larger or smaller than the partial 
derivatives suggest. In this section, we now consider the total 
effect on the calculated abundances, both from all the corre- 
lated changes in determining of the model atmosphere param- 
eters and also from the changes in the EWs of the elements 
being measured. 

For our particular linelist, T en -, log g and £ all affect the 
slope of the excitation potential vs. Fe I abundance plot, the 
slope of the reduced equivalent width vs. Fe I abundance plot 
and the difference in abundance between Fe I and Fe II lines. 
Therefore, using the partial derivatives we calculated above, 
we simultaneously solve for the changes in T e ff, log g and £ 
that resulted in the slopes of the two lines and the difference in 
the two ions being the same as in the unmagnified case. The 
changes in the atmosphere models from the original parame- 
ters in Table 1 are listed in Table 2. The changes are given 
as magnified parameters - unmagnified parameters. Finally, 
we interpolated a model with those parameters and run the 
EWs through that model to calculate the final change in the 
abundances. We checked that the final model has the excita- 
tion, equivalent width and ionization equilibrium to within the 
scatter (e.g., < 0.01 dex difference in ionization equilibrium) 
expected considering the effects of rounding off temperatures, 
gravities and velocities. As these are interpolated models to 
begin with, these changes will never be that precisely deter- 
mined. We note that the metallicity of the model atmosphere 
will have an effect on the abundances. However, the changes 
in loge(Fe) based on Fe I lines, which is how the metallicity 
of the atmosphere is usually set, are small enough that we will 
ignore this last effect for the purposes of the paper. For exam- 
ple, the change of 0.02 dex in Fe I for the 5600K model with 
Profile 35 would lead to a subsequent change for 0.006 dex in 
Fe I if folded back into a new model. Of course, the detailed 
analysis of a particular event could take these into account. 

Figure [lTI shows the resulting changes in the loge for the 



elements considered here for 5000K, 5600K and 6200K. This 
figure illustrates the small magnitude of the changes expected 
from DLM in the dwarfs. Also, the changes in abundances 
are not always in the same sense. Depending on what part of 
the source is most highly magnified, the change caused by not 
taking DLM into account can either be an increase or decrease 
over the true abundances in the star. For some ions, such as 
Fe II, the changes using the new parameters are smaller than 
using the new EW with the unmodified atmosphere (seen by 
comparing Figure 6 and Figure 11). For others, such as Fe I, 
the changes with the new parameters are larger than seen just 
using the new EWs. 

5. APPLICATION TO THE STUDY OF MICROLENSED 
BULGE DWARFS 

As stated in the Introduction, the emphasis on high magni- 
fication events means that in many cases there will be DLM 
when the spectra are taken. Our ability to determine the ef- 
fect on the spectra from observations falls into four differ- 
ent possibilities: 1) the lightcurve is indistinguishable from a 
point source and DLM is truly negligible; 2) the lightcurve 
is indistinguishable from a point, but DLM occurs; 3) the 
lightcurve shows extended source effects, which also means 
that DLM occurs; 4) there is no photometric data for crucial 
parts of the event. Of these possibilities, the first means that 
the spectrum can be analyzed using standard spectroscopic 
techniques, while the third can either be analyzed with stan- 
dard techniques, leading to errors of the order presented here, 
or analyzed using spectra with the DLM derived from the 
lightcurve. It is cases 2) and 4) that are of the most con- 
cern for ensuring reliable abundance results. Case 2) is fully 
addressed by Event A, which illustrates the maximum DLM 
that is possible without giving rise to noticeable effects, given 
a well-covered light curve. Case 4), i.e., poor light curve cov- 
erage, must be handled on a case-by-case basis. To date, all 
microlensing events with spectra have had adequate photo- 
metric coverage. 

6. CONCLUSION 

Microlensing of bulge dwarfs provides the exciting oppor- 
tunity to obtain high S/N and high-resolution spectra of stars 
that otherwise would be unfeasible with current telescopes. 
However, microlensing can produce finite source effects that 
affect the spectra by differentially magnifying annuli in the 
star compared to an unmagnified star. We have examined 
the practical impact of such DLM by synthesizing spectra 
of dwarfs with magnification profiles similar to events that 
have been observed and performing an abundance analysis. 
We find that, as expected, there are changes to the spec- 
trum, which results in changes in the atmospheric parameters, 
metallicities and abundance ratios. Given the small size of the 
changes, the S/N of the observed spectra and the accompa- 
nying error in the atmospheric parameters, and the length of 
the exposures compared to the duration of a specific magnifi- 
cation profile, the resulting effect on the abundances is small 
compared to other sources of error. We note that the change 
in abundances can be either positive or negative (see Profile 
35 vs. Profile 18) depending on the form of the magnification 
profile. Microlensing events are more commonly in the shape 
of Event A, with mild magnification causing limb brighten- 
ing, so there will be a tendency for events to be biased to 
lower metallicity,but Figure 1 1 shows that Event A leads to 
changes of < 0.01 dex in metallicity. These results, com- 
bined with the fact that finite source effects do not affect all 



9 



tie 
CO 

a 
3 



cm 
cd 



0.05 


0.05 
0.05 



0.05 
0.05 



0.05 



~i 1 1 r 



T 



~i 1 1 r 



~i 1 1 r 



~i 1 1 r 



T=5000K 



O 



o 



o 



□ n □ 



CI Ctl Til 
Sil Sell Til 



Oi00 0° 

A ▲ 

▲ A 



° o 
□ D 9i- 



o 



VI Crl Mnl Fell Nil Znl Zrl N d I h 
Crll Fel Col Cul Yll Cell 



H — I — I — \— t 
T=5600K 



H — I — I — h 



H — I — I — h 



H — I — I — h 



+ 



□ □ 



o 



°" UU Q n°flfi 00 n°fiO 
□ nHnOntfD □ nQo 



O 







▲ ▲ 



□ 



Q O- 



Cl Cal Til VI Crl Mnl Fell Nil Znl Zrl N d I H 
Sil Sell Till VII Crll Fel Col Cul Yll Cell 



H — I — I — h 



H — I — I — h 



H — I — I — h 



H — I — I — h 



+ 



T=6200K 



O n ^ Q o 0n OO Q OpoQ°Q 



□ 



□ 



□ 

▲ A 



□ 

o 



□ 



CI Cal Til VH 
Sil Sell Till 



Crl Mnl Fell Nil Znl Zrl N d I H 
Crll Fel Col Cul Yll Cell 



_1 I I L_ 



J_ 



J I I L_ 



J I I l_ 



_1 I I L_ 



J_ 



Element 



FIG. 1 1 . — The change in abundance (unmagnified-magnified) from the unmagnified case for three profiles studied in detail in this paper. Top: 5000K model 
Middle: 5600K case Bottom: 6200K case. Not all elements were measured for each temperature. We see that when DLM is not taken into account when 
analyzing the spectrum, the abundances can be in error by up to 0.05 dex. 



10 



TABLE 2 

Changes in Parameters for Magnified 
Model Atmospheres 

Profile AT cff Alog g A [Fe/H] A£ 



Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 

Profile 
Profile 
Profile 



35 



35 



35 



35 



35 



35 



35 



35 
Ori 



35 



Ori 



35 



35 



35 



Original T cff =5000K 

-4.0 -0.006 -0.005 0.002 

23.5 0.024 0.008 0.031 
-40.3 -0.039 -0.018 -0.059 

Original T cff =5100K 

-4.7 -0.006 -0.005 0.002 

23.8 0.028 0.008 0.036 

-31.9 -0.010 -0.015 -0.047 

Original T cff =5200K 

-5.6 -0.005 -0.001 -0.013 

21.7 0.025 0.007 0.035 
-39.6 -0.046 -0.008 -0.048 

Original T eff =5300K 

-4.7 -0.005 -0.001 -0.013 

20.6 0.025 0.008 0.032 
-40.2 -0.049 -0.018 -0.065 

Original T eff =5400K 

-4.5 -0.007 -0.003 -0.010 

20.8 0.022 0.005 0.041 
-36.9 -0.055 -0.017 -0.059 

Original T cff =5500K 

-4.4 -0.016 -0.001 -0.003 

19.8 0.020 0.007 0.039 
-34.5 -0.060 -0.016 -0.062 

Original T cff =5600K 

-4.2 -0.007 -0.004 -0.002 

21.2 0.026 0.008 0.041 

-35.7 -0.052 -0.020 -0.063 

Original T cff =5700K 

-5.0 -0.003 -0.003 -0.010 

20.1 0.026 0.006 0.045 

-38.2 -0.051 -0.017 -0.075 

Original T eff =5800K 

-6.8 0.015 -0.003 -0.018 

25.7 0.023 0.010 0.054 
-36.9 -0.047 -0.019 -0.079 

Original T eff =5900K 

-14.4 -0.043 -0.007 -0.013 

15.9 0.002 0.005 0.044 
-45.6 -0.094 -0.023 -0.073 
;inal T c ff=6000K, high gravity 

-4.0 -0.033 0.001 -0.011 

20.7 0.019 0.008 0.045 
-35.8 -0.059 -0.014 -0.071 
ginal T c ff=6000K, low gravity 

-6.9 -0.033 -0.002 -0.003 

18.8 0.026 0.004 0.042 

39.6 -0.063 -0.015 -0.072 
Original T cff =6100K 

-7.0 -0.010 -0.003 -0.010 

19.7 0.052 0.004 0.031 
40.1 -0.063 -0.019 -0.082 

Original T eff =6200K 

-4.9 -0.020 0.000 -0.014 

18.1 0.033 0.007 0.034 

39.2 -0.096 -0.019 -0.073 



of the high-magnification cases for which dwarfs have been 
observed, demonstrate that this is not the cause of any differ- 
ences between the observed metallicity distribution function 
for giants and dwarfs. Finally, if the lightcurve shows even 
more extreme finite source effects than have been modeled in 
this paper, magnification profiles for that event can be con- 
structed and model spectra calculated that appropriately take 
those profiles into account, so this will not be a limiting prob- 
lem for the accuracy of abundances derived for microlensed 
dwarfs. 

We thank Carlos Allende Prieto and Judy Cohen for useful 
comments and discussion. Thomas Masseron provided vital 
help in the setup and use of Turbospectrum. 
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